Accepted for publication in the Astrophysical Journal 

Preprint typeset using 1^1^^ style cmulatcapj v. 03/07/07 



X0-2b: TRANSITING HOT JUPITER IN A METAL-RICH COMMON PROPER MOTION BINARY 

Christopher J. Burke^ P. R. McCullough\ Jeff A. Valenti^, Christopher M. Johns-Krull^, Kenneth A. Janes^, 

J. N. Heasley", F. J. SuMMERs\ J. E. Stys\ R. Bissinger'"', Michael L. Fleenor", Cindy N. Foote'', 
Enrique Garci'a-Melendo**, Bruce L. Gary^, P. J. Howell^, F. Mallia"', G. Masi", B. Taylor^, T. Vanmunster^^ 

(Received 2007 Sep 14) 
Accepted for publication in the Astrophysical Journal 



o 
o 

(N 



6 



> 
cn 
o 
o 
p 

o 
i> 
o 



■ 

X 



ABSTRACT 

We report on a V=11.2 early K dwarf, XO-2 (GSC 03413-00005), that hosts a i?p=0.98±H? Rj, 
Mp=0.57± 0.06 Mj transiting extrasolar planet, X0-2b, with an orbital period of 2.615857±0.000005 
days. XO-2 has high metallicity, [Fe/H]=0.45±0.02, high proper motion, fitot = 157 mas yr~^, and has 
a common proper motion stellar companion with 31" separation. The two stars are nearly identical 
twins, with very similar spectra and apparent magnitudes. Due to the high metallicity, these early 
K dwarf stars have a mass and radius close to solar, M^, = 0.98 ± 0.02 Mq and i?* = 0.97±q o^ Rq. 
The high proper motion of XO-2 results from an eccentric orbit (Galactic pericenter, Rper < 4 kpc) 
well confined to the Galactic disk {Z^ax ~ 100 pc). In addition, the phase space position of XO-2 
is near the Hercules dynamical stream, which points to an origin of XO-2 in the metal-rich, inner 
Thin Disk and subsequent dynamical scattering into the solar neighborhood. We describe an efficient 
Markov Chain Monte Carlo algorithm for calculating the Bayesian posterior probability of the system 
parameters from a transit light curve. 

Subject headings: binaries: eclipsing - planetary systems - stars: individual (GSC 34130-0005) — 
techniques: photometric - techniques: radial velocities 



1. INTRODUCTION 

We announce the discovery of an extrasolar planet, 
X0-2b, that transits a bright, V=11.2, star. With an 
orbital period, P ~ 2.6 days, planetary radius, Rp =0.98 
Rj, and planetary mass, Mp =0.57 Mj, X0-2b belongs 
to t he growing class of transi ting Hot Jupiter (HJ) plan- 
ets (jCharbonneau et ani2007l) . Despite the increasing 
number of transiting planets known, much work still re- 
mains to understand the observed properties of p lanets. 
The transi ting planets HD209458b (Charbonn cau et al.l 
[2000; Hcn rv et al.ll2000D. HAT-P-lb (Bakos c t aI1[200l. 
and \Y ASP-lb (jCameron et al.ir2007l : ICharbonneau et alj 
l2007b( l have anomalously large radii compared to theo- 
retical models and are thought to require a n external 
source of energy to rema in inflated (Bodcnh eimer et al.l 
l2003t iGuillot fc Showma n 2002: Win n fc Holmanll2005D . 
However, recently iBurrows et al., (2003) explain the radii 
of transiting planets without the need for an extra source 
of energy by accounting for enhanced metallicity opaci- 
ties and properly comparing observed radii to theoreti- 
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cal radii. In contrast, HD149026b has an extremely high 
density and small radius (iSato et al. 2005) . There is gen- 
eral agreement that a pure H/He mixture cannot explain 
the small radius of HD 149026b, but a planet model with 
a massive central core, Mc ~ 70 M®, of heavy elements 
along with a small H/He envelope ca n explain the ra dius 
of HD149026b (Sato et a l. 2005; For tnev et al.ll2006[ ). 

Disentangling the effects of stellar irradiation, mi- 
gration, central core mass, and composition on the 
observed properties of plan ets requires discovering 
more bright transiting planets (jCharbonneau et al.ll2007t 
iFortnev et~al. 2006). X0-2b is the second contribution 
to the brigh t transiting planet samp le provided by the 
XO Project (iMcCulloug h et aLl[2005l): XO-lb being the 
first (|McCullough e t al. 200l). Of the bright transiting 
planet hosts {V < 12), XO-2 has the highest metallic- 
ity ([Fe/H]=0.45 see § 13. ip . In addition, from a recent 
catalog of extrasolar planets ([Butler et al.ll2006l ). XO-2 
has a higher metallicity than 96% of all known extrasolar 
planet hosts. 

Metallicity plays a crucial role in the formation of plan- 
ets and the resulting HJ atmospheres. The frequency 
of radial velocity detected planets is kno wn to increase 
with metallicity ([Fischer fc Valenti[ |2005( ). In the core 
accretion planet formation model, a high metallicity en- 
vironment grows larger cores and enables more objects 
to reach the critical mass necessary for runaway gas ac- 
cretion and transforma t ion into a detectable gas giant 
planet ([Ida fc Linll200l iBenz et all [20061) . Detailed fits 
to the mass and radius of the known transiting planets 
yield s a mass estimate for the central refracto ry element 
core ([Guillot et al.ll2006t[&j^rows et al.ll2006f ). These in- 
vestigations derive a larger core mass for planets that 
transit higher metallicity stars. In addition to the bulk 
properties of extrasolar planets, metallicity plays an im- 
portant role in the planet's atmosphere, especially for a 
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HJ experiencing large stellar irradiation. An increased 
metallicity results in a greater absorption of the stel- 
lar irradiation and larger equilibrium tem peratures than 
a com parable planet of lower metallicity (jFortnev et al.l 
[2OOI . The variations in equilibrium temperature can 
lead to variations in the dominant observable features in 
the planet's atmosphere. 

Amongst the other transiting HJ planets, X0-2b 
shares a commo n characteristic w ith another transiting 
HJ, HAT-P-lb (jBakos et al.ll2007D . Both planets orbit 
one member of a nearly-equal-mass, wide separation stel- 
lar binary. The XO-2 stellar binary system has a sepa- 
ration of 31" 4600 AU separation with a distance 
of ~ 150 pc) and both components have an identical 
within the uncertainties 157 ma s yr~^ prop er motion vec- 
tor as measured with Tycho-2 (iHog et al.i,2000) a nd the 
high proper motion catalog of ijLepine fc Sharal I2005D . 
It is not unusua l for p lanets to exist in binary systems. 
iRaghavan et al.l (|2006[ ) find >23% of stars with planets 
have a stellar companion, however, they find evidence 
that the extrasolar planet sample is deficient in stellar bi- 
naries when compared to the field. A stellar companion 
with a 4600 AU separation is not expected to influence 
the planet formation process. Even assuming an orbital 
eccentricity, e=0.8, for the unknown eccentricity of the 
XO-2 stellar binary, pl anets within Oc ~ 170 AU of XO-2 
are dynamically stable (jHolman fc Wieger3ll999t ) (we as- 
sume the current measured projected binary separation 
is close to the actual semi- major axis of the orbi t when 
calculating Oc). Empirically. iDesidera fc Barbieril (|2007l ) 
do not find any statistically significant difference between 
planets around single stars and planets in wide (a > 100 
AU) separation binaries. However, understanding the 
influence of a stellar companion on the planet formation 
process is complicated by the fact that the binary con- 
figuration during planet formation may be vastly differ- 
ent than what is currently o bserved (jMarzari fc Barbieril 
[2007t iMalmberg et al.ll200l . 

The similar brightness of the XO-2 stellar binary com- 
ponents and their angular separation provides an ex- 
cellent opportunity for detailed line abundance stud- 
ies. The higher metallicity of the known planet host 
stars is thought to be of primordial origin rather 
than heavy element po l lution due to i nfalling planets 
(|Fischer fc Valentil [2005t FGonzalezI [200I . The infalling 
planet model for explaining the higher metallicity of the 
known planet host stars has predominately come from 
observations showing differing element abundances be- 
tween common proper motion binary components. These 
element abundance differences have typ ically been over- 
turned by further independent analyses ([Gonz alez 2006) . 
The difficulty in reliable abundance differences results 
from systematic effects in analyzing stars of different T^f / 
and evolutionary state (jSchuler et al.ll2006h . Such effects 
are reduced for the XO-2 stellar binary. 

On the sky, XO-2 and a line joining its stellar compan- 
ion has a position angle of 162°, nearly along the North- 
South direction (see the finder chart in Figured]). XO-2 is 
the Northern component of the binary (indicated by the 
arrow in Figure [1]) . The magnitude difference between 
the components is AV — 0.05 mag. Thus, in the optical, 
XO-2 is fainter than its Southern companion. The tra- 
ditional nomenclature for binary stars designates XO-2 
as X0-2B and the planet X0-2Bb. However, throughout 
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Fig. 1. — Finder chart for XO-2 and its stellar companion at 
two epochs. The arrow indicates XO-2. North is toward the top of 
the figure and East is toward the left of the figure. (Left) Digital 
Sky Survey POSS I Red image from 1953. (Right) 2MASS H-band 
image from 2000. The common proper motion is predominately 
directed South. 

this article we simply designate the transiting HJ stellar 
host XO-2 and the planet around it X0-2b. When it is 
unclear from the context we alternatively designate XO-2 
as X0-2N and the Southern stellar companion as X0-2S. 
We adopt this naming convention because the magnitude 
difference between the stellar components is small and 
the spectra are nearly identical, both of which make it 
difficult to distinguish the objects, whereas the 31" sepa- 
ration readily distinguishes the stellar components based 
on declination. We refer to both stellar components as a 
unit as the XO-2 stellar binary. 

In § [21 we provide details of the discovery and follow 
up photometry and spectroscopy. The high resolution 
spectroscopy is analyzed in § 13.11 to determine the stel- 
lar properties of X0-2N and X0-2S. The high precision 
photometry is analyzed in § 13.21 employing an efficient 
Markov Chain Monte Carlo (MCMC) algorithm to de- 
termine the properties of X0-2b. We confirm the plan- 
etary mass of X0-2b with radial velocity measurements 
in § 13.31 The ephemeris for X0-2b is refined and the 
transit observations are investigated for transit timing 
variations in § 13.41 There is a brief discussion regarding 
the XO-2 stellar binary and its Galactic orbit in §[51 and 
we summarize the article in § [51 

2. OBSERVATIONS 

2.1. XO Project Photometry 

X0-2b is the second transiting planet discovered by 
the XO survey af te r XO- lb ([McCullough ct al. 2006). 
'McC ullough"etaI] ((20051 ) and iMcCuUough fc Burki 
(2007) describe the instrumentation, operation, and 
analysis in more detail than the summary provided here. 
The twin, 200mm XO cameras power drift scan in dec- 
lination over a 7 deg x 62 deg strip every ten minutes on 
clear nights for more than 2 months per season of visi- 
bility. The XO observations employ a broad (0.4/im to 
0.7/im) passband. The star XO-2 comes from the XO 
strip centered on RA 8.0 hr and is one of several thou- 
sands of bright (V< 12) stars monitored in this strip 
over the period September 2003 to September 200 5. The 
Box- fitting Least Squares algorithm, BLS, (Kovacs "et al.l 
[2OOI was employed to search the nearly 3000 observa- 
tions per star for repetitive transit events with periods 
ranging from P = 0.5 — 10 day. We perform the transit 
search on two realizations of the light curve. One re- 
alization is the ca l ibra ted light curves as described in 
iMcCullough et aT (2005) and the other realization of 
the li ght curve has the SysRcm algorithm ([Tamuz et al.l 
l2005f ) applied to the calibrated light curve in order to fur- 
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Fig. 2. — Discovery phased light curve from the XO Project data. 
(Top) The phased light curve around the transit event. {Bottom) 
The phased light curve over the full orbital period. The individ- 
ual measurements (points) are shown binned (solid line) to reduce 
noise. The transit occurs at At = 0.0 day and any secondary 
eclipse for a circular orbit would occur at At = 1.3 day (dotted 
line). No secondary eclipse is evident above the noise. At the res- 
olution of the XO cameras, both components of the XO-2 stellar 
binary are within the photometry aperture causing dilution of the 
transit depth. 

ther remove systematics. In the case of XO-2, there were 
no substantive improvements in the hght curve quality 
after implementing SysRem and both analyses identified 
the same transit period and ph ase. 

McCuUough fc Burk^ ()2007[ ) describe the selection of 
transit candidates for followup. In addition, for consid- 
eration th e candidates n iust p ass selection criteria as de- 
scribed in iBurke et al.l (|2006D . In brief, we require ob- 
servations covering more than 1.5 transit events, avoid 
0.5 and 1.0 day periods where systematic aliases result 
in false-positive detections, require the transit depth, 
Sm < 0.1 mag, and require the transit to have higher 
significance than systematic errors (as mea sured by the 
transit to anti-transit Ratio Statistic, RS, of lBurke et al.l 

dlooe)). 

Figure [2] shows the XO light curve phased on the de- 
tected period as returned from BLS. We achieved 0.7% 
or 0.007 mag per observation precision on this V=11.2 
object. The transit occurs at At = 0.0 day and any sig- 
nature of a secondary eclipse would occur at At — 1.3 
day for a circular orbit. At the angular resolution of the 
XO survey, both components of the XO-2 stellar binary 
are within the photometry aperture leading to transit 
depth dilution. Light curves resolving the components of 
the XO-2 stellar binary result in twice the transit depth 
(see Figure [3]), but in the discovery light curve the tran- 
sit depth is 0.007 mag. With the XO cameras, we ob- 
served two nearly complete transits of X0-2b at the end 
of 2004, and beginning of 2005 (Julian dates 2453355 and 
2453376). In addition, on seven other nights, XO cam- 
eras captured partial ingress and egress events. First 
occurrence of a partial transit was at the end of 2003 
on Julian date 24552994. Table [1] provides a sample of 
the photometry for XO-2 from the XO cameras. The full 
table is available in the online edition. 

2.2. Extended Team and Follow Up Photometry 

The Extended Team (E.T.) provides photometric fol- 
low up for XO candidates. The E.T. (R. B., M. F., C. F., 
E. G-M.,B. G., P. H., F. M., G. M., and T. V.) is a 
collaboration of professional and amateur astronomers 



(jMcCullough et al.ll2005l I2006D . We sent the candidate 
list containing XO-2 to the E.T. on January 16, 2007. 
On, January 18, 2007, E.T. observations confirmed the 
transit events for X0-2b in the R-band using a 0.6-m 
telescope (solid green line in Figure [3|) . The initial ob- 
servations were truncated due to high air mass, but the 
observations confirmed the XO ephemeris, the transit oc- 
curs in the Northern star of the XO-2 stellar binary, and 
the undiluted transit depth is consistent with a planetary 
radius and the XO photometry. A complete transit was 
observed by the E.T. on January 26, 2007 and the binned 
R-band hght curve obtained with a 0.35-m telescope is 
shown as the solid red line in Figure [3l The light curves 
shown in Figure [3] are used for refining the ephemeris 
and looking for evidence of transit timing variations in 
§ 13.41 Based on observations of other transit candidates 
by the E.T., repeat transit events in the same passband 
have a 0.2% or 0.002 mag standard deviation in deriving 
the transit depth. At this level of precision, the E.T. 
observations of X0-2b are consistent with a gray tran- 
sit. Table [1] provides E.T. photometry for XO-2. For the 
E.T. light curves, the median differential magnitude out 
of transit provides the flux normalization and the stan- 
dard deviation out of transit provides the uncertainty in 
the measurements. 

On February 16, 2007 we observed a transit event of 
X0-2b with the 1.83-m Perkins Telescope at Lowell Ob- 
servatory using th e PRISM instrument in imaging mode 
panes et al.ll2004l ). The PRISM camera is a 2048x2048 
Fairchild CCD with 0.39"/pix resolution. The transit 
event was well positioned in the evening, occurring over 
the airmass range 1.05< X <1.16. To improve efficiency 
a subframe containing X0-2N and X0-2S was read out 
with a 8-10 s cadence. The average seeing during the 
course of observations, 3", was poor for the site. The final 
light curve is a differential light curve for X0-2N using 
X0-2S as the comparison shown in Figure |4] employing 
aperture photometry and an R-band filter. The selected 
aperture size for photometry minimizes the resulting rms 
scatter in the differential light curve. In addition to the 
R-band photometry, photometry in the BVI passbands 
was obtained before the transit, at mid-transit, and after 
the transit. The gap in the R-band data at mid-transit 
accommodates the multiple filter data. A computer fail- 
ure resulted in the gap of R-band data at the start of 
ingress. 

To normalize the differential light curve, the average 
magnitude from out of transit data was subtracted from 
the light curve. This was done independently for each 
half of the light curve separated by the mid-transit gap. 
There is a Am = 0.001 magnitude difference in the nor- 
malization zeropoints between the data before and after 
the transit. This offset results from repositioning the 
telescope and refocusing for the BVI data acquisition at 
mid-transit. During an uninterrupted R-band data se- 
ries, a 0.3 pix rms positional accuracy was maintained 
resulting in ^^0.0014 mag rms over 1 minute intervals. 
Accurate positioning was not maintained for the BVI 
data resulting in ~ AO. 005 mag offsets between the data 
obtained before and after the transit. The average BVI 
passband data before transit was subtracted from the 
BVI data obtained at mid-transit and after transit to 
yield the differential light curve in these passbands as 
shown in FigurelH Within the systematics resulting from 
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Fig. 3. — Binned light curves from the Extended Team for XO- 
2. The color of the text in the lower left corner indicates the date 
(2007), observer, and passband of the observations. The passband 
labels are Johnson for the system except Clr indicates unfiltered 
observations and BB indicates a blue blocking (> 0.5 /im) filter. 

repositioning of the stars on the detector, the transit is 
gray. 

We obtained photometric B, V, Rc, and Ic magni- 
tudes for X0-2N and X0-2S using a 0.35-m telescope 
on the photometric nig ht of January 2 4, 2007 (Tabled]). 
A single Landolt area (|LandoltllT99l was observed at 
the same airmass as XO-2 to derive the zero point and 
color transformation coefficients. The color transfor- 
mation coefficients were consistent with comprehensive 
standard star measurements from four Landolt fields ob- 
tained two weeks previous using the identical instrumen- 
tal setup. The color range of Landolt standards was 
-0.14 < (B - V) < 1.4. The B, V, Rc, and Ic absolute 
photometric accuracies are 0.04, 0.04, 0.04, and 0.05 mag 
r.m.s., including both the rms scatter around the photo- 
metric transformation model and an estimated system- 
atic error. The Tycho-2 magnit udes for XO- 2 listed in 
Table [2] transform (via Table 2 of lBesseU (j2000D ) to John- 
son V = 11.25, i.e. 0.07 mag (2-ct) fainter than our es- 
timate. In addition, we accurately measured the instru- 
mental magnitude difference between X0-2N and XO- 
2S. We find X0-2N is fainter than X0-2S by 0.07±0.008 
mag in the B-band, 0.055±0.004 mag in the V-band, 
0.040±0.004 mag in the R-band, and 0.030±0.003 mag 
in the I-band. 

2.3. Spectroscopy 

After confirmation of the XO transit light curve 
from E.T. observations (see W2.2\\ . we initiated queue 
schedule observations of X0-2N and X0-2S with the 
High-Resolution Spectrograph (H RS), a fiber fed cross- 
dispersed echelle spectrograph ()TulH |T998), on the 
McDonald Observatory 11-m Hobby-Eberly Telescope 
(HET) in order to measure the mass of the planet. The 
first HRS observations using an iodine gas cell for pre- 
cision radial velocities commenced on January 26, 2007. 
Table [3] provides dates of the HRS observations along 
with the resulting radial velocities. The instrument setup 
provides R=60,000 resolution and wavelength coverage 
over the range 4000 < A < 7800 A with center at A = 
5900 A. We extracted the two-dimension al echel l e spec - 
tra using procedures described in Hinkl e et al.l (|200ClD . 
The resulting Signal-to-Noise ratio (SNR) varied from 
20-50 per extracted pixel at the blaze peak. We calcu- 
late radial velocities for X0-2N and X0-2S in § [SJl 
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Fig. 4. — Light curve from the 1.8m Perkins Telescope at Low- 
ell Observatory in the R-band (points) for XO-2N using XO-2S 
as the comparison star. The open symbols show the differential 
photometry in the B (diamond), V (triangle), and I (square) pass- 
bands using the average of the BVI observations obtained before 
the transit for the flux normalization level (dotted line). Accurate 
positioning of the stellar image on the detector was not maintained 
during the BVI observations resulting in 0.5% or 0.005 mag sys- 
tematic offsets. However, accurate positioning of the stellar image 
on the detector was maintained during the R-band light curve re- 
sulting in ~0.0014 mag rms over 1 minute intervals. Also shown is 
the best-fit transit model in a sense during the MCMC analysis. 
(solid line). 

To measure the stellar parameters of X0-2N and XO- 
2S, we also o btained spectra with the 2dCoude echelle 
spectrometer (i Tull et al."l995!) on the McDonald Obser- 
vatory 2.7-m Harlan J. Smith Telescope (HJS). We ob- 
tained two spectra of X0-2N and a single spectrum of 
X0-2S with R=60,000 and wavelength coverage of 3900- 
9600 A. We determine the stellar parameters in § 13.11 

3. ANALYSIS 
3.1. Stellar Properties 

We use the Spectroscopy Ma de Ea sy (SME) analysis 
pack age of Valenti & Piskunovl (|1996f ) with refinements 
from iValenti fc Fischer. (.20051 ) on the HJS spectra to 
measure the stellar properties of X0-2N and X0-2S. 
We briefly describe the process here. The free param- 
eters, Teff, logg, [M/H], vsini, [Na/H], [Si/H], [Ti/H], 
[Fe/H], and [Ni/H], are varied in order to minimize the 
difference of the resulting synthetic spectrum to the ob- 
served spectrum. A quadratic continuum is fit over 8 
wavelength intervals for each unique set of the above 
free parameters. When generating the synthetic spec- 
trum, the pressure-temperature profile of the atm osphere 
comes from interpolating the atmosphere grid of IKurucd 
(|1992D . The atomic fine hst comes from the Vienna 
Atomic Line Database (VALD) (iPiskunov et all Il995f) 
and the molecular line list comes from iKuruc 
The line strengths and van der Waals damping parameter 
of the line list database were adjusted to improve agree- 
ment with the ob s erved solar spectrum, as described in 
IValenti fc Fische^ (l2005h . 

Table|3]lists the stellar parameters for X0-2N and XO- 
2S based on the SME analysis. The results for X0-2N 
are based on the average from two spectra and the results 
for X0-2S are based on a single spectrum. Both XO- 
2N and X0-2S independently result in a metal enhanced 
abundance, [Fe/H]=0.45±0.02. Using the primary ob- 
servables from the SME analysis (Te//, abundances, and 
logg) and the apparent magnitude in the V-band ('ii l2.2p . 
we determine secondary stellar properties, M^, R^, dis- 
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tance, and age using the isochrones (lYi et al. l2001h 
following the procedure of lValenti fc Fisched (|2005l ) . The 
distance to XO-2 is unknown, thus the probability den- 
sity function for M*, i?^,, and age are calculated for a 
sequence of trial distances in steps of 10 pc. The SME 
isochrone analysis for select distances to XO-2 are listed 
in Tabled We determine a distance, d=150 pc, to XO-2 
in § 13.21 from a joint analysis of the SME analysis and 
transit light curve. We show the probability density for 
the X0-2N parameters in Figure O for the preferred dis- 
tance to XO-2 of 150 pc. 

The spectrum alone yields an estimate of the stel- 
lar gravity, loggsme- Estimates of Mi, and Ri, from 
the isochrone analysis provide an additional estimate of 
\oggiso as a function of distance to XO-2. The condition 
loggsme—^oggiso yiclds an approximate distance to XO- 
2. The data for X0-2S provides a consistency check, and 
ideally the condition loggsme =log(7iso for X0-2S is met 
for the same distance as X0-2N. In practice, we find XO- 
2S has loggsme = 4.6; too large a value for loggiso to ac- 
commodate at any distance. Also, X0-2S being brighter 
and thus more massive than X0-2N implies loggsme for 
X0-2S should be lower than loggsme = 4.5 of X0-2N, op- 
posite of what is measured. Analysis of the transit light 
curve along with the physical parameters for X0-2N from 
the isochrone analysis in §[321 yields loggidso ~ 4.5, con- 
firming the spectroscopic logf^s^e for X0-2N is correct 
and the loggsme for X0-2S is an overestimate. 

The overestimate in loggsme for X0-2S is not without 
precedent. The SME based logg can result in higher val- 
ues than the isochrones allow for some objects in the 
Spectroscopic Properties of Cool Stars (SPOCS) cata- 
log as a result of the numerous degeneracies that ex- 
ist between the s p ectra l parameters (see Figure 16 of 
IValenti fc Fischeil (|2005f )). Despite individual cases of 
an overestimated logg, comparing stars analyzed with 
SME in common with other independent spectral analy- 
ses d oes not reveal a sys t emati c offset in logg (see Figure 
20 of lValenti fc Fischeil (|2005[ )). The 1-ct uncertainties 
in the stellar parameters given in Table |4] are based on 
the typical rms scatter in parameters measured in inde- 
pendent, multiple spectra for stars in the SPOCS cat- 
alog. Based on the 1-cr uncertainties in Table IH the 
difference in logg between X0-2N and X0-2S is signifi- 
cant. However, the distribution of repeat measurements 
for parameters in the S POCS catalog (see Table 5 of 
IValenti fc FischeJ ()2005D ) display extended wings and 
the Gaussian-based uncertainties underestimate the pos- 
sibility of outlying measurements. Thus, we also pro- 
vide the 99.7% con fidence intervals f o r the parameters 
based on Table 5 of lValenti fc Fisched (|2005D . Since, the 
X0-2S parameters are based on a single spectrum, the 
99.7% confidence interval is more appropriate than rely- 
ing solely on the 1-cr error when deciding on the signifi- 
cance of any differences between the properties of X0-2N 
and X0-2S. 

3.2. Markov Chain Monte Carlo Light Curve Analysis 

iFordI ()2005f ). iGregorvl ()2005f ). and references therein 
provide a thorough discussion of the theory behind 
Markov Chain Monte Carlo (MCMC) Bayesian analysis 
along with a practical MCMC implementation for radial 
velocity planet detection. In a multidimensional prob- 
lem, the MCMC algorithm is an efficient means of calcu- 
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Fig. 5. — Distributions for four stellar parameters derived from 
the SME analysis (see text) for a distance of 150 pc. The values 
of the mean and limits containing ±34% of the distribution from 
the mean are annotated on the figures and listed in Table |4] along 
with corresponding values for distances of 140 pc and 170 pc. 



lating the Bayesian p osterior probability for parameters. 
iHolman et al.l (|2006l ) describe MCMC analysis applied to 
determining the system properties for transiting extra- 
solar plane ts. In their an alysis of the transit light curve 
for XO-lb, Ho lman et al.l (|2006ll calculate the posterior 
probability using the stellar radius, i?^, planet radius, 
Rp, orbital inclination, i, and the two coefficients of the 
quadratic limb darkening law as free parameters. In the 
transit fitting problem, the relationship between R^,, Rp, 
and i has a high degree of degeneracy and nonlinearity. 
We illustrate the degeneracy and nonlinearity between 
parameters in the top panels of Figure |S| for a MCMC 
calculation of the X0-2b system properties. The nonlin- 
ear, "banana-shaped" degeneracy between Rp and i slows 
the rate of convergence for MCMC algorithms. MCMC 
algorithms work more efficiently when the relationship 
between parameters is multi-normal without covariance 
(jKosowskv et al.ll2002D . 

To remove the nonlinearity between parameters, we 
choose the following set of free parameters: i?^,, p = 
Rp/Ri,, and the total transit duration from I''' to 4* 
contact, T. With this nonlinear transformation, the de- 
generacy between the new set of parameters is along a 
straight line as shown in the middle panels of Figure [HI A 
further linear transformation between parameters yields 
an eigenbasis set of parameters with a multi-normal non- 
covariant relationship that results in an efficient MCMC 
calculation. These transformations improve the rate of 
convergence for the MCMC analysis by more than a fac- 
tor of 100 when measured by the autocorrelation length 
of the samples in the MCMC sequence. Appendix [X] 
provides further details of the MCMC implementation 
developed for this study. 

We employ the Metropolis-Hastings algorithm with a 
normal proposal distribution for calculating the Markov 
Chain. We follow a Gibbs-like sampling technique where 
each step in the chain consists of a number of intra steps 
updating each individual parameter in turn. Several 
short, trial chains iteratively yield scale factors of the 
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Fig. 6. — Correlation between parameters for samples in a 
MCMC calculation of the XO-2b posterior distribution. The pa- 
rameters ij*, Rp, and i are tightly correlated and the Rp versus i 
correlation is nonlinear (Top). With the choice of parameters, i?*, 
p, and T, the correlations between parameters are linear (Middle). 
Solving for the eigenvectors of the covariance matrix from samples 
shown in the Middle Panel yields an eigenbasis set of parameters 
(Lower). The eigenvector along the direction of largest variance, 
ei is predominately along ij*. Eigenvectors 62 and are a mix- 
ture of p and T. The multi- normal without covariance properties of 
the eigenbasis set of parameters provides a more efficient MCMC 
calculation. 



normal proposal distribution for each parameter with a 
25% to 40% acceptance rate for the trial samples. The 
likelihood function is given by e~^'^^ , were we have as- 
sumed the errors are normally distributed, and the data 
have uniform weights, is the squared difference be- 
tween observations an d the analytic transit model of 
iMandel fc Agol (I200I ). The model assumes negligible 
eccentricity. The observations are the R-band data from 
the 1.8m Perkins Telescope shown in Figure [H 

The calculation has seven free parameters: Mi,, Ri,, p, 
T, to, ui, and U2. to is the mid-transit time offset from the 
ephemeris with a period given by the XO observations 
and a mid-transit zeropoint near the mid-transit gap of 
the light curve, HJD 2454147.75. The hmb darkening 
coefficients, ui and U2, model the limb darkening with 
the quadratic law, / = 1— ui(l— /x)— M2(1— /^)^, where / is 
the specific intensity normalized to unity at the center of 
the stellar disk and /x is the cosine of the angle between 
the line o f sight and the surf ace normal. In practice, 
we follow iHolman et "aD ()2006f ) by adopting ai = ui -I- 
2u2 and 02 = 2ui — U2 as the parameters used in the 
calculation. This linear combination of limb darkening 
coefhcients reduces their mutual degeneracy. 

The prior for i?* , p, and r is given by Equation IA5[ 
which is equivalent to a prior uniform in Ri,, Rp, i. The 
priors for Ri,, p, and to have cutoff values well beyond 
values allowed by the data. The prior on r has an upper 
limit cutoff, T < Tmax, where Tmax is the longest transit 
duration possible for a given Ri,, M^,, and p. We assume 
uniform priors for the limb darkening coefficients with 
the following physically motivated limits on the param- 
eters. We require the highest surface brightness to be 
located at the disk center (ui > 0.0), require the specific 
intensity to remain above zero (ui + U2 < 1-0), and do 
not allow limb-brightened profiles (ui -I- 2u2 > 0.0). The 
form of the prior for is a Gaussian where the cen- 
tral value and standard deviation of the Gaussian are a 
function of Ri, as given by the SME isochrone analysis 
data given in Table [5] and described in § 13.11 We employ 
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Fig. 7. — Marginalized posterior probability for the XO-2 and 
XO-2b parameters from MCMC samples. We adopt the median 
of the posterior probability (arrow) as the point estimate of the 
parameter. The solid lines indicate the 68.3% credible interval. 

a spline interpolation over the grid of Mi, and gm as a 
function of i?^. 

The estimate of the posterior probability comes from 
7 independent chains of length Nchn = 60000 with vary- 
ing initial conditions. This results in an effective length 
Neff = Nchn/Ncor = 12000 after taking into account the 
autocorrelation length Ncor = 5. Using these 7 chains, 
the largest Gelman-Rubin statistic amongst the param- 
eters, R— 1.0 002, where R<1.02 ind icates convergence 
of the chain ijGelman fc RubinI [l99l ) . Figure [7] shows 
the resulting posterior probability distribution for each 
parameter after marginalization over the other parame- 
ters. The posterior probability is simply a normalized 
histogram of the MGMG sample values. We adopt the 
median as the best single point estimate of the poste- 
rior probability. To derive an a% credible interval for 
a parameter, the N MCMC samples are sorted by the 
parameter of interest. The lower limit of the credible 
interval is taken as the ((1 — a)/2)N^^ sorted sample, 
and the upper limit of the credible interval is taken as 
the (1 — (1 — a)/2)N^^ sorted sample. The arrow point 
along the abscissa in Figure [7] indicates the median of 
the posterior probability, and the the vertical solid lines 
in Figure [7] show the 68.3% credible interval. 

Figure [8] shows the joint posterior probability for the 
stellar limb darkening coefficients, ui and U2- The solid 
contours are isoprobability contours containing 68% and 
90% of the MCMC samples. The smaU points illustrate 
the remaining 10% of the samples lying outside the re- 
gion of highest probability. The region of highest proba- 
bility for the limb darkening coefficients differs from the 
theoretically calculated limb darkening coefficients (open 
triangle symbol) obtained from Claret (2000). The dash 
lines illustrate the prior limits for the limb darkening co- 
efficients. The lower dash line in Figure [5] corresponds to 
not allowing limb brightened specific intensity profiles. 

Using the SME isochrone analysis, we translate the 
MCMC samples for i?^ into distance to XO-2 and age 
estimates. Similar to the procedure that defines the 
prior on M^,, we interpolate the R^ versus trial distance 
and age relationships given by the SME isochrone anal- 
ysis (i.e., for a given stellar radius estimate, the SME 
isochrone analysis has a best distance and age estimate 
for the system). The posterior probability for these two 
parameters is shown in Figure O We repeated the cal- 
culation assuming a prior uniform in cosi, and this did 
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Fig. 8. — Joint posterior probability for the stellar limb darken- 
ing coefficients, ui and U2. The solid contours are isoprobability 
contours containing 68% and 90% of the MCMC samples. The 
remaining 10% of the samples lying outside the region of highest 
probability are also shown (points). The region of highest proba- 
bility differs from the theoretically calculated limb darkening coef- 
ficients (triangle). The prior limits for the limb darkening coeffi- 
cients are indicated with dashed lines. 

not materially affect the parameter estimates. We sum- 
marize the properties of X0-2b in Table [T] 

The uncertainties given in Table [7] for X0-2b and Ta- 
ble[2]for XO-2 represent the precision of our experiment. 
These uncertainties represent the expected scatter of val- 
ues obtained if the experiment was repeated with similar 
quality data and identical procedures. Other systematic 
sources of error affect the accuracy of our measurement 
that are most likely comparable or larger than our pre- 
cision. The sources of systematic error only enter into 
our prior for M^,. Our adopted uncertainty in AI^, ^^2% 
follows from the uncertainty in T^ff and metallicity in 
the SME isochrone analysis (see § I3.ip aided by the very 
we ak dependence of on the unknown distance to XO- 
2. iHillenbrand fc Whitd (|2004D show the isochrones 
employed in this study successfully predict stellar masses 
to within l%-3% for main sequence stars Af^ > 0.6 
Mq with independent dynamical mass estimates. For 
the derived Teff a i id [Fe /H] for XO-2, the isochrones 
from iGirardi et all ()2OO20 at maximu m differ by 3% in 
Mi, from the isochrone prediction. ICodv fc Sasselovl 
(11)02) find 7% systematic error in for HD 209458 
due to uncertainty in Helium abundance and the treat- 
ment of convection. In light of these potential sources 
of systematic uncertainty, we increased the standard de- 
viation of the Gaussian prior on to cr = 0.07 Mq. 
This larger uncertainty on M^, resulted in more symmet- 
ric and slightly broader posterior distributions for Ri, 
(cr±0.03 Rq), Rp (cr± 0.03 Rj), and age (cr± 1.4 Gyr). 
The surface gravity of the p lanet, g^, is independent of 
Mi- ([Southworth et all |2007( ) . From the measured ra- 
dial velocity semi-amplitude and light curve parameters, 
gp =14.8±1.6 m s~^ for X0-2b, where the uncertainty in 
gp is dominated by the uncertainty in the radial velocity 
semi-amplitude. Other parameters weakly dependent on 
Mi, and available from the light curve observations di- 
rectly are a/R^, = 8.2±'^-l and a/Rp = 79.0±0;9. 

3.3. Radial Velocity Measurements 

We measured the mass of X0-2b usi ng the radial veloc- 
ity technique described below and bv 'McCulloug h et al.l 
(|2006) . The 2-dimensional spectra obtained with the 
HET were extracted to 1-dimensional spectra and as- 
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Fig. 9. — The marginalized posterior probability for the XO-2 
that rely on the SME isochrone analysis. The SME isochrone anal- 
ysis translates the posterior distribution for ij* into the posterior 
distribution for Distance and Age. 

sociated approximate wavelength solutions derived from 
ThAr spectra obtained during twilight (§ 12. 3p . To ac- 
count for optical distortions, the wavelength solution fits 
the centroids of thousands of ThAr lines to a function of 
the X and Y coordinates of the CCD that includes terms 
linear in X and in Y plus the following cross terms: XY, 
XXY, XYY, XXYY, XXXY, and XYYY. Using a down- 
hill simplex minimization algorithm, "Amoeba," we 
adjusted the parameters of a synthetic spectrum to fit, in 
separate ~15 A sections, the stellar spectrum observed 
through an iodine absorption cell. The requirement of 
strong iodine lines limits radial velocity estimates to the 
wavelength range 5210 < A < 5700 A. The synthetic 
spectrum consists of a high-re solution spectrum o f the 
Sun, the Earth's atmosphere^'^ (Wall ace et al.lfT998h. and 
a hig h-resolution spectrum of an iodine gas cell ( CochranI 
|2000| ) convolved with a Voigt profile to approximate the 
line-spread- function of the instrument. In addition to 
the convolution, we used a few additional free parame- 
ters to model specific physical or instrumental character- 
istics: the radial velocity of the star, a shift of the iodine 
spectrum attributable to instrumental errors in the ap- 
proximate wavelength solution specific to the particular 
15 A section, the continuum level, and an exponent that 
adjusts the optical depths of the solar absorption lines in 
order to better approximate the (non-solar) stellar spec- 
trum. The difference of the parameters representing the 
radial velocity of the star and that attributable to the 
instrument from the iodine spectrum equals the topocen- 
tric radial velocity of the star, which we transform to the 
barycentric frame of the solar system. We average the 
individual stellar radial velocity estimates from each of 
the 15 A sections to determine the stellar radial veloc- 
ity measurement at each epoch and its associated 1-a 
uncertainty (see Table [3]). 

Figure [10] shows the resulting radial velocity curve 
phased with the X0-2b ephemeris determined from the 
transits and assuming zero eccentricity. The typical un- 
certainty for each measurement, any = 20 m s^^. The 
radial velocity semi-amplitude, K=85±8 m s~^. This 
amplitude resuffs in Mp =0.57±0.06 Mj for X0-2b, as- 
suming M^=0.98 Mq for XO-2 and a circular orbit for 
X0-2b. The six radial velocities measured for X0-2S 

^•^ Within the wavelength range of interest, telluric absorption 
lines are negligible but are included anyway. 
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Fig. 10. — a) The radial velocity of XO-2 oscillates sinusoidally 
with a semi-amplitude K = 85it8 m , implying XO-2b's mass is 
0.57it0.06 Mj. b) The period and phase of the radial velocities were 
fixed at values determined by the transits. The mean stellar radial 
velocity with respect to the solar system's barycenter has been 
subtracted. In order to determine K, we used the EIET spectra 
calibrated with an iodine absorption cell (filled circles), c) In this 
representation of the data, a circular orbit yields a straight line of 
slope — K. 

(Table [3]) are consistent within the observational errors 
with no radial velocity variation and show no significant 
evidence for a HJ orbiting X0-2S. 

3.4. Ephemeris and Transit Timing Variations 

The database of E.T. transit light curves for X0-2b 
enable us to refine the ephemeris for X0-2b and to quan- 
tify transit timing variations from this ephemeris. The 
transit light from the 1.8m Perkins Telescope on Febru- 
ary 16, 2007 provides a precise measurement of the mid- 
transit zeropoint of X0-2b. From the best fitting transit 
model in a sense during the MCMC calculation, we 
find the ephemeris zeropoint is 2454147.74902 with an 
uncertainty of 17 s. 

To refine the orbital period of X0-2b, we calculate mid- 
transit times for the nearly complete transit events ob- 
served with the XO cameras on 2453355 and 2453376 
HJD and the E.T. transit events shown in Figure [3] with 
nearly complete ingress and egress coverage. We com- 
pare the transit data with a limb-darkened transit model. 
The only free parameter of the transit model is the mid- 
transit time. The other parameters of the transit model 
are fixed at the best fitting transit model in a sense 
during the MCMC calculation, and the model employs 
theoretically calcul a ted li mb darkening coefficients ob- 
tained from IClaretl ()2000f ) in the appropriate passband 
of the observations. The transit model has a flux decre- 
ment of p^/2 when applied to the XO survey data to 
account for the transit depth dilution caused by the flux 
of X0-2S in the photometric aperture. 

Table [S] provides the resulting mid-transit times. The 
three independent mid-transit estimates on 2454168 HJD 
enables estimating the typical uncertainty (l-cr=3 min) 
in the transit timing. To refine the orbital period we 
minimize the difference between the transit tim- 
ing observations and ephemeris model. We keep the 
ephemeris zeropoint fixed at the value derived from the 
1.8m Perkins Telescope light curve. The best-fit period, 
2. 615857±0. 000005 day, results in a xL« = 10.1 for 
1^=11 degrees of freedom. This suggests there are no 
significant timing variations. For the minimization we 
assumed a uniform error, a— 3 min, on the transit tim- 



ings. The ephemeris for X0-2b accumulates a 5 min 
uncertainty by 2010. 

4. DISCUSSION 



The transit candidate of'Mandushev et al.' (2005') illus- 
trates the non-negligible potential for triple stars to have 
transit light curves and radial velocity variations consis- 
tent with a planet. However, in the case of X0-2b, our 
attempts to explain the light curve and spectroscopy with 
a physical stellar triple fail. We employ the isochrone 
appropriate for the the physical properties of XO-2 sup- 
plemented with the low-r nass stellar isochrone between 
0.072< <0.5 M0 from lChabrier et al.l (|2000D . stellar 
limb darkening coefficients from IClaretl i200W and the 
light curve synthesis routine of IWilsonr()1993[ ) to model 
a stellar triple system. The constraints on the transit 
duration and transit depth from the light curve require 
> 0.95 Mq for the primary in a stellar binary blended 
with the light of XO-2. The required stellar binary has 
>75% the fiux of XO-2 and has a radial velocity semi- 
amplitude, K > 16 km s~^ . Such a binary would be read- 
ily apparent in the spectrum of XO-2 given the narrow 
spectral features for XO-2, v sinz < 2.3 km s~^. We can- 
not completely rule out the possibility of a line-of-sight 
faint background binary blended with the light of XO- 
2 as an explanation for the observations. However, the 
sinusoidal shape of the radial velocity variations necessi- 
tate the line-of-sight binary to have a systemic velocity 
similar to XO-2 otherwise the radial veloc ity curve de- 
velo ps asymmetries that are not observed ([Torres et al.l 
|2005 f). Given the large proper motion of XO-2, ground- 
based adaptive optics or space-based observations have 
the potential to definitively rule out a background line- 
of-sight binary. 

In addition to XO-2, only HAT-P-1, HD 20782, HD 
80606, HD 99492 (GJ 429B), HD 178911B, and HD 
186427 (16 Cyg B) have <2.0 m ag difference between 
the stellar binary components (|Desidera fc Barbieril 
120071 supplemented by data from the SIMBAD 
database) and can be considered nearly-equal-mass stel- 
lar binaries hosting a known extrasolar planet. In 5 out 
of the 7 binary systems, the planet orbits the lower mass 
star of the binary. In the remaining 2 systems, HD 20782 
and HD 80606, the planet orbits the more massive star, 
but the orbits for these planets have the highest eccen- 
tricities amo ngst all known plan ets, e=0.93 and e=0.92, 
respectively ( Butler et al.ll20()a) . Equivalently, if we re- 
strict the sample of extrasolar planets in nearly-equal- 
mass stellar binaries to modest eccentricity (e<0.8), in 
5 of the 5 known systems, the planet orbits the lower 
mass star. We offer several speculative possibilities to 
explain the current pattern of planets in nearly-equal- 
mass stellar binaries: 1) The pattern arises due to small 
sample statistics, 2) Planets form more readily around 
the lower mass component of a nearly-equal-mass stellar 
binary, 3) Planets form equally likely around the compo- 
nents of a nearly-equal-mass stellar binary, but planets 
orbiting the primary component have lower detectability 
due to lower pla net mass or higher orbital eccentricity 
(|Cumming|l200l . 

The XO-2 stellar binary has a higher space ve- 
locity with respect to the Local Standard of Rest 
(LSR), V = VU^~+WTW^ - 100 kms~ \ than most 
of the other known extrasolar planets (jSantos et al.l 
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120031 lEcuvillon efall l2007t ) . According t o the p urely 
kinematic cl a ssifica tion of iBensbv et al.l ()2003f ) and 
iBensbv etall (|2006[ ). XO-2 has ~7 times higher proba- 
bility of belonging to the Thick Disk than the Thin Disk. 
From a chemical abundance perspective, a metallicity of 
[Fe/H]< —0. 3 is typically the upper limit for Thick Dis k 
metallicit ies (iMishenina et alj 120041 : iReddv et all I2006D . 
However. IBensbv et al.l ( 2006() has recently claimed kine- 
matically selected Thick Disk stars with [Fe/H]~ 0.0 
have a chemical abundance pattern distinguishable from 
Thin Disk stars at the same metallicity. The sample of 
stars with kinematics suggestive of the Thick Disk but 
Thin Disk metallicities, like XO-2, is still too small to de- 
termine their relationship to th e tradition al Thin / Thick 
Disk populations. However, 'Misheni na et al.l ()2004f ) 
distinguishes between objects with Thick Disk kine- 
matics and Thin Disk metallicity by determining the 
maximum height above the Galactic disk their orbit 
attain s, Zmax- Both the estimate of Bens bv et al.l 
(|2005f ) for Zmax and integrating the orbit of XO-2 us- 
ing the axi-symm e tric, s tatic Galaxy potential model of 

1 Allen fc SantillanI (|199H ) yield Z^ax ~ 100 pc for the or- 
bit of XO-2. Given the much larger scale height for the 
Thick Disk, aTP ^ 100 kpc, under the classification of 
IMishenina et all (|2004D . XO-2 belongs to the Thin Disk, 
since its orbit is well confined to the Galactic plane. The 
high space velocity of XO-2 results from high eccentric- 
ity rather than exc ursions from Galac t ic pla ne. Using the 
simple model from I Allen k, SantillanI (|1991f ). the orbit of 
XO-2 has a pericenter within ~4 kpc of the Galactic cen- 
ter and apocenter ^9 kpc (e~0.4). 

Currently, XO-2 lags the LSR by F = -78 kms'^ and 
is moving away from the Galactic center at J7 = —71 
kms^^. This places XO-2 in pha s e space near the Her- 
cules stream (jFamaev et al.ll2005t lEcuvillon et al.ll2007[) 
{U = 42, au = 28, y = -52, cry = 9 kms"!), but is 
lagging the LSR more than the typical Hercules stream 
member. The possibility that XO-2 belongs to a dynam- 
ical stream originating from the metal-ric h inner Galaxy 
is an example that supports the findings of lFamaev et al.l 
(|lb05) and Ecuvillo net al.l ([2007) that compared to field 
stars, metal-enhanced dynamical streams should be over- 
abundant in detectable extrasolar planets. 

5. SUMMARY 

The star XO-2, GSC 34130-0005, hosts an approxi- 
mately Jupiter-size, 0.6 Jupiter-mass transiting extra- 
solar planet, X0-2b, with an orbital period ^2.6 days. 
XO-2 is a V=11.2 early K dwarf with high metallicity, 
[Fe/H]=0.45, high proper motion, jitot = 157 mas yr~^, 
and has a common proper motion stellar companion 
with 31" separation. The foUowup high resolution spec- 
troscopy yields M^=0.98±0.02 Mq for the mass of XO- 

2 and Mp=0.57±0.06 Mj for the mass of X0-2b. The 
foUowup high precision photometry yields _R^=0.97±q q^ 
Rq for the radius of XO-2 and Rp^QM±l f^ Rj for 
the radius of X0-2b. Joint analysis of the light curve 
and spectroscopy yields an isochrone based age of XO- 
2 t=5.3±J y Gyr, and an isochrone based distance to 
XO-2 (i=150±2 pc. The quoted values and their un- 
certainties are Bayesian credible intervals encompassing 
68% of the marginalized posterior probability. For the 
Bayesian analysis we assume a Gaussian-like prior on 



Mi,, derived from the spectroscopic and isochrone anal- 
ysis of XO-2 and uniform priors on i?^, i?p, orbital in- 
clination, and limb darkening coefficients. We describe 
an efficient Markov Chain Monte Carlo method to cal- 
culate the Bayesian posterior probability for the system 
parameters from a transit light curve. 

X0-2b adds to the sample of Jupiter-mass planets re- 
siding in nearly-equal-mass stellar binaries. In 5 of the 
7 nearly-equal- mass stellar binaries hosting a planet, the 
planet orbits the lower mass stellar component. In the 
remaining 2 stellar binaries where the planet orbits the 
higher mass star, the planets' orbits are highly eccentric 
(e>0.9). Equivalently, if we restrict the sample of extra- 
solar planets with modest eccentricity (e<0.8) orbits in 
nearly-equal-mass stellar binaries, in 5 of the 5 known 
systems, the planet orbits the lower mass star. We spec- 
ulate on possible astrophysical reasons for this pattern 
beyond the simple fact that the statistics are based upon 
a very small sample and thus may be spurious. 

With its high proper motion, XO-2 has kinematics sug- 
gestive of Thick Disk membership but a Thin Disk metal- 
licity. In contrast to the typical Thick Disk member, the 
high proper motion of XO-2 results from an eccentric or- 
bit (Galactic pericenter, Rper < 4 kpc) well confined to 
the Galactic disk, Z„iax ^ 100 pc. XO-2 may originate in 
the metal-rich inner Thin Disk and was dynamically scat- 
te red into the sol ar neighbor hood. Similar to t he fin dings 
of lFamaev et al.l (|2b05. ) and lEcuvillon et al.l (|2007D . the 
discovery of X0-2b suggests metal-enhanced dynamical 
streams from the inner Galaxy may be abundant in de- 
tectable extrasolar planets. 
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APPENDIX 

MARKOV CHAIN MONTE CARLO DETAILS 

Ca re must be taken to properly assign prior probabilities for a nonlinear transformation in parameters ()Chu et al.l 
l2003f ). In particular, a uniform prior in i is not the same as a uniform prior in t. The transformation law of 
probabilities provides the necessary prior probability to maintain the uniform prior in i?^-, i?p, and i for our chosen set 
of free parameters: i?*, p, and r, where 



P . R. 
— arcsm — 
TT \ a 



(1 + - (a/i?^ cos i)^ 



1 



cos 



1/2N 



(Al) 



where P is the orbital period, a is the semi-major axis, and we have assumed zero eccentricity (jSeager fc Mallen-Ornelai 
[20031) . 

The general transformation law with multiple dimensions states 



p(yi,?/2, ...)dyidy2-- ^p{xi,X2, ■■■) 



d{xi,X2, ...) 



d{yi.V2, ■■■) 



dyidy2-., 



(A2) 



where the original joint probability distribution, p{xitX2t ■■■)dxidx2---, is transformed into another probability dis- 
tribution in terms of the new set of variables (yi,2/2, ■••) by multiplication with the absolute value of the Jacobian 
determin ant, T he new y variables must have the same number and be expressible in terms of the old x 

variables ([Press et al.lll993 ). Writing the old variables in terms of the new variables, R^, ^ R'^, Rp = pR'^, and 



arccos < 



'{R'Jani + pY 



smc 



(1 



dR 
dR 



1 



dRp 
dR' 
di 
dRL 



OR*. - 
dp - 
dR^ 



di 







dR^ 
dr 
OR,, 



where = ttt/P, the Jacobian matrix is written 

\\dO/dO\\ = 

where we have indicated the trivial elements with value one or zero. Overall, the Jacobian simplifies to 



= 
= 



di 



(A3) 



(A4) 



\{dRp/dp){dl/dT)\ 



-R' sin ( 



{{RUay{i + py 



1) 



(sin02-(i?;/a)2(l + p)2)- 



(A5) 



With our choice of parameters, the proper prior to maintain uniform priors in Ri,, Rp, and i is the above Jacobian 
multiplied by the original joint probability distribution for uniform priors, which is a constant. The prior probability 
goes to zero for < r < Tmax, where Tmax = P/tt arcsin (i?*/a)(l -I- p). 

A trial chain is necessary to define the covariance matrix for the final eigenbasis set of parameters. Mi,, i?*, p, r, 
ui, and U2 enter into the covariance matrix. We find to shows no significant correlation with respect to the other 
parameters. To determin e the necessary hnear transformation, we determine the eigenvectors of the covariance matrix 
built from a trial chain ([Tegmark et al.|[2004D . This is equivalent to a Principal Component Analysis. The bottom 
panels in Figure [6] show trial samples in terms of the eigenvectors of the covariance matrix. The eigenvector along the 
largest variance, ei, is predominately along Ri,, and 62 and 63 are a mixture of p and r. 

The autocorrelation length of samples in the MCMC sequence is one method to quantify the efficiency of the 
calculation. With the original set of parameters (i?*, Rp, and i), the correlation length (when the autocorrelation 
drops by half) in the i?* parameter varied, 500 < Ncor < 1500 steps. Using our new set of parameters (i?*, p, and r), 
40 < Ncor < 60 steps. Finally, with the eigenbasis parameters, 4 < Ncor < 5 steps, where the correlation is measured 
in the physical variable i?*. 
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TABLE 1 
XO Survey & E.T. Light Curve Data'' 



Heliocentric Julian Date 


Light Curve 
[mag] 


Uncertainty 
(l-cr) [mag] 


Filter 




Observatory 


2452961.14380 


-0.0027 


0.0039 


W 


1 


XO 


2452961.14404 


-0.0002 


0.0037 


W 


1 


XO 


2452964.11621 


-0.0057 


0.0034 


W 


1 


XO 


2452964.11646 


-0.0018 


0.0033 


W 


1 


XO 


2452964.12329 


0.0034 


0.0034 


W 


1 


XO 



^ The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 
^ Average of N measurements 



TABLE 2 
Stellar Properties 



Parameter 


XO-2N 


XO-2S 


Reference 


CSC ID 


03413-00005 


03413-00210 


a 


RA (J2000.0) 


7'"48™06^.47 


7*" 48^07^. 48 


a,b 


Dec (J2000.0) 


+50°13'33".0 


-l-50°13'03".3 


a,b 


Galactic Latitude b [dcg] 


29.33 




a 


" Longitude 1 [dcg] 


168.29 




a 


V 


11.18±0.03 


11.12±0.03 


c 


(B-V) 


0.82±0.05 


0.79±0.05 


c 


(V-Rc) 


0.49±0.05 


0.46±0.05 


c 


(V-Ic) 


0.86±0.05 


0.82±0.05 


c 


Vt 


11.24 


11.20 


b,d 


(B-V)t 


0.70 


0.86 


b,d 


J 


9.74±0.02 


9.74±0.02 


e 


(J-H) 


0.40±0.03 


0.37±0.03 


e 


(H-K) 


0.03±0.03 


0.10±0.03 


e 


Spectral Type 


KOV 


KOV 


c 


Distance [pc] 


150±4 




c 


Ha [mas yr~^] 


-34.7±2.6 


-33.1±2.9 


b 


^lS [mas yr^i] 


-153.6±2.4 


-154.1±2.7 


b 


Total /x[mas yr"-*^] 


157 


158 


f 


Radial Velocity (Bary) [km s~^\ 


47.4±0.5 




c 


U Space Velocity [km s~^] 


-72.0 




b,c,g,h 


V " 


-78.0 




b,c,g 


W " 


-4.6 




b,c,g 


Stellar Mass [Mq\ 


0.98±0.02 




c 


Stellar Radius [Rq\ 


n 07+0.02 




c 



References: 

a) SIMBAD 

b) Tycho-2 (H0e et al. 2000) 

c) this work 

d) Tycho-2 photometry on the Johnson system ^Bessell 2000; ^amaje^_et_^l^ 2002) 

e) 2MASS (Skrutskie et al. 2006) 

f) LSPM (Lepine fc Shar a 2(M^ 

g) Space Velocity w.r.t. LSR after correction for Solar peculiar motion (Mihalas_&_^inne^^98^. 

h) Negative U is away from the Galactic Center 
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XO-2N 



XO-2S 



TABLE 3 
Radial Velocity Shifts 



Object Julian Date Radial Velocity U ncertainty 

-11 



-245000 


Shift [m s-i] 


(l-cr) [m s 


4127.6626 


-60.6 


22 


4128.6459 


99.5 


18 


4134.8318 


-34.0 


15 


4135.8366 


-57.8 


19 


4136.6494 


71.0 


19 


4136.8477 


64.6 


20 


4158.7782 


-78.5 


19 


4159.7882 


73.6 


20 


4160.5891 


27.3 


23 


4160.7715 


37.0 


26 


4133.6383 


-21.6 


18 


4134.6257 


-1.5 


15 


4135.8481 


-32.4 


20 


4136.8599 


-7.6 


17 


4158.7932 


30.6 


17 


4168.7572 


16.7 


17 



TABLE 4 
Results of the SME Analysis 



Parameter 


XO-2N 


XO-2S 


Uncertainty (l-cr) 


Uncertainty (99.7%) 




5340 


5500 


32 


233 


logg [cm s ] 


4.48 


4.62 


0.05 


0.36 


11 sin i [km s^-"-] 


1.4 


1.2 


0.3 


2.1 


[M/H] 


0.44 


0.45 


0.02 


0.20 


[Na/H] 


0.49 


0.63 


0.02 


0.18 


[Si/H] 


0.39 


0.47 


0.02 


0.12 


[Ti/H] 


0.36 


0.42 


0.04 


0.26 


[Fe/H] 


0.45 


0.47 


0.02 


0.22 


[Ni/H] 


0.50 


0.52 


0.02 


0.16 


[Si/Fe] 


-0.06 


0.00 


0.03 


0.25 



TABLE 5 

Spectroscopically Derived Stellar parameters 



Parameter @ 140 pc @ 150 pc @ 170 pc 





0.98 


0.97 


0.96 


Mass [Mq] 


0.99 


0.98 


0.98 


1.00 


1.00 


0.99 




0.90 


0.95 


1.07 


Radius [Rq] 


0.91 


0.97 


1.10 


0.93 


0.98 


1.12 




4.49 


4.42 


4.32 


Log(g) [cm s-2] 


4.51 


4.45 


4.34 


4.53 


4.48 


4.37 




0.77 


3.55 


8.35 


Age [Gyr] 


2.04 


5.27 


9.49 




3.68 


6.92 


10.62 



For each parameter, the middle row is the maximum likelihood value, and the 
values in the rows above and below span the 68% likelihood of the probability 
distributions (cf. Figure [SJl. The three columns correspond to three assumed 
distances for XO-2. 
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TABLE 6 
Mid-Transit Times 



Heliocentric Julian Date'' Observatory ID' 



-2450000 


3355.14209 


XO 


3376.07349 


xo 


4126.82324 


MF 


4126.82324 


BG 


4134.66992 


MF 


4142.51660 


EM 


4147.75049 


BG 


4160.82812 


BG 


4168.67578 


BG 


4168.68066 


OF 


4168.67871 


MF 



^ Based upon independent observations on 
2454168 HJD, uncertainty in mid-transit 
time, a=3 min. 

^ Observatory ID is author initials, except 
XO is XO survey data. 



TABLE 7 
The Planet XO-2b 



Parameter Value 

P 2.615857±0.000005 d 

tc 2454147.74902±0.0002 (HJD) 

K 85±8 m s-i 

a 0.0369±0.0002 A.U. 

i 88.9±0.7 deg 

Mp 0.57±0.06 Mj 

Rp 0.98±g »f Rj 

gp 14.8±1.6m s-^ 



